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Abstract 

We analytically derive the effective layer elastic energy of smectic-A and smectic-C liquid crystals 
by adiabatic elimination of the orientational degree of freedom from the generalized Chen-Lubensky 
model. In the smectic-A phase, the effective layer bending elastic modulus is calculated as a 
function of the wavelength of the layer undulation mode. It turns out that an unlocking of the 
layer normal and the director reduces the layer bending rigidity for wavelengths smaller than 
the director penetration length. In the achiral smectic-C phase, an anisotropic bending elasticity 
appears due to the coupling between the layer displacement and director. The effective layer 
bending rigidity is calculated as a function of the angle •& between the layer undulation wave- 
vector and the director field. We compute the free energy minimizer •d = 9. It turns out that 
6 varies from 0° to 90° depending on the tilt angle, undulation wave-length and other elastic 
constants. We also discover a new important characteristic length and the discontinuous change 
of 9. Using the elastic constants of Chen-Lubensky model, we determine the parameters of the 
more macroscopic model [Y. Hatwalne and T. C. Lubensky, Phys. Rev. E 52, 6240 (1995)]. We 
then discuss the hydrodynamics, and demonstrate the alignment of director and the propagation 
of the anisotropic layer displacement wave in the presence of an oscillatory wall and a vibrating 
cylindrical source respectively. 
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I. INTRODUCTION 



Liquid crystals have many fascinating physical properties, because they have both the 
molecular position and the molecular orientation as degrees of freedom [1]. This multiplicity 
gives rise to interesting phase transition phenomena, including the tricritical and the Lifshitz 
point in the phase diagram [2], where the first and second order transitions are switched [3, 
4], and three phases (i.e., the nematic, smectic-A (Sm-A) and smectic-C (Sm-C)) coexist [5, 
6], respectively. To explain the phase behavior around the nematic-Sm-A-Sm-C (NAC) 
Lifshitz point theoretically, a variety of models have been developed [7-12]. De Gennes [7] 
and Chu and McMillan [8] introduced a tilt angle and an in-plane director as new degrees 
of freedom, respectively. Chen and Lubensky [9] took fourth order derivative of the density 
field into account. Benguigui [10] and Huang and Lien [11] adopted a smectic-C scalar order 
parameter together with a smectic-A order parameter. Grinstein and Toner [12] also used 
the two-dimensional in-plane director component to describe the free energy. Among them, 
the Chen-Lubensky model is supported by quite a few X-ray scattering experiments [13-15]. 

In the smectic phase, the layer order and the director often cause a frustration. In chiral 
liquid crystals, geometrical incompatibility of a uniform smectic order and a helical director 
configuration leads to a rich variety of defect phases [16]. They can be discussed with the 
Landau-de Gennes model, utilizing an analogy with superconductors in the magnetic field, 
where the smectic order parameter, the director vector and the chirality are identified with 
the wave function of a superconducting particle, the electromagnetic vector potential and 
the external magnetic field, respectively [17]. The twist-grain-boundary (TGB) phase is 
the simplest defect phase where the groups of planes, containing parallel screw dislocations, 
are regularly stacked with the dislocations in adjacent planes tilted each other at a constant 
angle, and finite length smectic slabs are inserted between the dislocation planes [18, 19]. The 
Chen-Lubensky model, a higher order extension of the Landau-de Gennes model, is again 
successful in explaining the structure and phase transition of the TGBa and TGBc phases, 
which have the Sm-A and Sm-C slabs respectively [20, 21]. More complex defect structures 
such as cholesteric and smectic blue phases have been discovered [16]. The cubic smectic 
blue phase (Sm-BPi) is comprised of a three-dimensional cubic disclination lattice, whose 
lattice constant is in the order of the wavelength of the visible light [22, 23]. The isotropic 
smectic blue phase (Sm-BPi so ) has fluid-like rheological properties and the detailed structure 
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is completely unknown [24]. Such novel defect structures in principle, should be understood 
with the Landau-de Gennes and Chen-Lubensky models presented above. However, in fact, 
these phenomenological models are not appropriate for theoretical explanations of these 
phases, because of the highly complex spatial structure together with the quite intricate 
free energy model. DiDonna and Kamien make use of a simpler phenomenological model to 
discuss the stability of the cubic blue phase [23]. In their Landau- Peierls form, the director 
is assumed to be parallel to the layer normal. However, a tilt of the director from the layer 
normal is reported to weaken the layer bending elasticity in the numerical simulation of 
the TGBa phase [25] , and to change the structural symmetry of the smectic blue phases in 
the experimental study [26]. Such anomalous properties in the TGB and the smectic blue 
phases have their origin in a nontrivial elastic mechanism of the undulated smectic layer, 
combined with the director degree of freedom. Therefore, a modification of the Landau- 
Peierls model, accommodating the director elasticity of Sm-A and Sm-C phases, would help 
the understanding of the complex layered structures. In addition, the Landau-Peierls model 
has the same layer compression and bending elastic energy forms as those of the other layered 
materials such as block copolymers and surfactant system [27, 28]. Thus, the characteristic 
features of liquid crystals are not very apparent in the original Landau-Peierls model. The 
director in liquid crystals should affect the layer elasticity, while block copolymers do not 
possess an orientational degree of freedom especially in the weak segregation limit [29]. Thus 
it is worthwhile to derive an intermediately simple model directly from the Chen-Lubensky 
free energy, eliminating the variation of the director degree of freedom, to clarify the role 
of the director elasticity in smectic liquid crystals. At the same time, by doing so, we can 
obtain the phenomenological elastic constants of the macroscopic models [23, 30], in terms 
of the more microscopic parameters of Chen-Lubensky model. 

Dynamics of Sm-C layers is also an interesting topic of liquid crystals [1, 31-34]. The 
director component parallel to layer (c-vector) plays an important role on the static and 
dynamic pattern formation due to the anisotropic elasticity and interaction with flow 
field [35, 36]. Shear flow orients the c-vector and results in a novel target pattern with 
a disclination [36, 37]. On the other hand, because of the static coupling between the layer 
displacement and the director, together with the rotational viscosity, an oscillatory wave 
traveling perpendicularly to the layer normal also rotates and aligns the c-vector in a cer- 
tain direction, even in a uniform oscillatory wave [1, 38]. It might be useful for an application, 
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for instance, mechanically-active optical devices and sensitive acoustic sensors [38-41]. 

This paper is organized as follows. In Section II, we review the original Landau-de Gennes 
model and derive the effective elastic energy of Sm-A phase. The analysis is extended to 
Sm-A and Sm-C phases with the generalized Chen-Lubensky model in Section III. Next we 
conduct a hydrodynamic simulation for Sm-C layer in Section IV. We conclude in Section 
V. 



II. LANDAU-DE GENNES MODEL FOR SMECTIC-A PHASE 

In this section, we briefly review the Landau-de Gennes model for Sm-A phase, and then 
calculate the effective layer elastic energy in terms of the layer displacement field. The 
physical meaning of the results are discussed qualitatively. 

Smectic liquid crystals can be described by the density modulation ^(r) and the director 
n(r). The density field \l/(r) is a complex order parameter with the absolute value being 
the amplitude of the smectic order and the phase describing the layer displacement [17]. 
The analogy between liquid crystals and superconductors leads to the phenomenological 
Landau-de Gennes model for the Sm-A liquid crystals 

F = -Flayer + -^cpl + ^Prank, (la) 



-^layer / dv 



'rm 2 + -\^\ A 



F cpl = J dr 1 - [B\\ \(n ■ V - iq ) + B ± \n x W| 2 ] , 



(lb) 



(lc) 



Fprank = J dr^ [K^V ■ Tlf + K 2 (n • V X nf 

+K 3 (n x V x nf] . (Id) 

The double-well potential Fi aycr describes the Nematic-Sm-A (NA) transition. F cp i and F-pr^ 
are the coupling energy between \1> and n, and the Frank elastic energy. The dimensionless 
temperature r is positive (negative) above (below) the NA point. The constants B\\ and B^ 
are the layer compression elastic coefficients which adjust the layer width to the equilibrium 
value d = 2n/qo. The Frank elastic constants -fQ (i — 1, 2, 3) correspond to the splay (i = 1), 
the twist {i = 2) and the bend (i = 3) deformations, respectively. 

We consider a perturbation of the uniform Sm-A structure far below the NA point to see 
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the effect of the director tilt from the layer normal. For this purpose, we set B\\ = B± = B 
and K\ — K 2 — K 3 = K to simplify the discussion. A more general case is considered in 
the next section for the Chen-Lubensky model. The amplitude of the smectic order is close 
to the equilibrium value at low temperature, so we assume 



tt(r) = exp[iq (z - u{r))} = 



(2) 



where the z-axis is set along the equilibrium layer normal, and the layer displacement field 
u(r) is introduced. The director n can be divided into its spatial average and the deviation: 

n(r) = e z + Sn(r). (3) 

The free energy components F cp i and Fp ra nk are now simplified as 

K f . ,_ „ |2 



Fcpl = 2A2 J dr \ Vu + 5n \ > ( 4a ) 
^ank = f / dr\VSn\ 2 , (4b) 

where A = K / B / q ) is the penetration length [18], and |V5n| 2 = (diSnj)(di5nj) 
(repeated indices are summed up). We next write the linearized effective free energy in 
terms of u(r). To do this, the director is adiabatically eliminated with the equation, 

5F 

(1 - nn) ■ -f- = 0. (5) 
on 

The factor 1 — nn ensures the normalization |n| 2 = 1, and 5n z = because \e z + 5n\ = 1. 
Using (4a), (4b) and (5), we obtain the Fourier representations 

5nj(q) = -k(q)iqju(q). (6) 



Thus the effective free energy is 

F _ K 
F cS - - 



u{q)\\ (7) 



where we define the z-axis and in-plane (x- and y-) axes as the parallel and perpendicular 
direction respectively (Fig.l). Equation (7) is identical with the Landau- Peierls free energy 
except for the Lorentzian dependence of the dimensionless layer bending elastic coefficient 
k(q), and the higher order cross term [2, 39]. 



The g-dependence of k(q) is interpreted qualitatively as follows. We consider a pure un- 
dulation case with a constant layer thickness q\\ =0. This condition is, in fact, a sufficiently 
good approximation to describe experiments in the thermodynamic limit [39]. Thus q \\ Sn 
is satisfied, which means that only the splay part \q ■ 5n\ 2 is nonzero in the three Frank 
elastic contributions. Since the Landau-de Gennes model (la) has no layer bending term in 
the ^/-dependent part, the effective elasticity comes solely from the splay term. In fact, if 
the director n coincides with the layer normal m, the splay term is turned into the layer 
bending elastic energy: 



K 
~2 



J dr{V ■ nf = 2K J drH 2 . (8) 
Here we denote the mean curvature of the smectic layer by H. 



(a) q ± « A" 1 (b) q ± > A" 1 

*n<-» !!!!!!!!!!!!!!!! III! 



I ( x, y, I HUM WM MM I ml I!! MM I (t! 

FIG. 1: Schematic representation of the undulated layer structure of Sm-A phase. The director 
and layer normal are (a) almost parallel for q± <^ A -1 and (b) decoupled for q± > A -1 . 

Equation (8) is achieved well for g^ 1 ^> A, where the director is locked to the layer normal 
(Fig. 1(a)). Oppositely, for gj 1 < A, the director cannot follow the layer deformation, and 
the splay term less contributes to the layer elasticity (Fig. 1(b)). We note that k(q) can be 
written also as a function of the root mean squared (RMS) tilt angle a = ^/ (a 2 ). With a 
smooth undulation u{r) = uq cos (q ± ■ r) (q± A -1 ), 

«(«) = 7273- (9) 

1 + (V2a\/u ) 2/3 

The detailed derivation is given in Appendix A. The effective layer bending rigidity decays 
with a. It implies that the coupling between m and n is certainly essential. 
In fact, for any the locking term is found in the coupling term (lc) 
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F cp i = - I dr 



n • -iW + ^o™(|V</>| - g ) + ^o|V0|(ra - n) 



n x ^-iW + ^o|V0|m^ 



(10) 



where the complex order parameter is decomposed into the amplitude and the phase: \&(r) = 
^o(r) exp(i(j)(r)) , and the gradient of the phase component defines the layer normal vector 
m(r) = V0/|V0|. In the parallel part, the first term favors a uniform smectic modulation 
amplitude, the second is the layer compression energy and the last term locks the director 
along the layer normal. In the perpendicular part, on the other hand, the first term is 
again the amplitude homogenizing contribution, and the second term reduces the director 
component perpendicular to the layer normal (c-vector) as long as the coefficient B± is 
positive. The Sm-C phase is stable with a negative B±, as we shall see in the next section. 



III. GENERALIZED CHEN-LUBENSKY MODEL FOR SMECTIC-A AND -C 
PHASES 

In this section, we introduce the Chen-Lubensky model in the generalized way. Then 
the effective elastic energy is calculated for both Sm-A and Sm-C phases. After a physical 
interpretation of the Sm-A energy, the effective layer bending elasticity of the Sm-C phase 
is discussed. The anisotropy is characterized by the angle d between the c-vector and the 
layer bending direction. The state diagrams for the easiest bending angle 6 are studied. The 
free energy as a function of d is also calculated. Finally, the model parameters of the more 
macroscopic free energy [30] are determined with those of Chen-Lubensky model. 



A. Generalized Chen-Lubensky model 

By adding higher order gradient terms, one can extend the Landau-de Gennes free energy 
to reproduce the achiral Sm-C phase. The Chen-Lubensky model is given in the most 
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generalized form: 



F 
' F 



Frank 



cpl 



cpl 



Frank 



2 r r 



r|^| 2 + -|^F| 
2 



C\\nirij + 



D|| InjUjA-Di^l 



JliTlj 



DiDjV) (SkDkD&Y + c.c. 



y rfr [(l + «i)(V-n) 2 
+ (1 + K 2 )is: 2 (n- V x n) 2 
+ {n x (V x n)} 2 ] , 



(11a) 
(lib) 



(11c) 



(lid) 



where Di = di—iqoUi is the covariant derivative and 5^ = Sij—riiUj is the projection operator. 
The elastic constants are C\\ and C± for the second order gradient term F^j, D\\, D\\± and 
D± for the fourth order terms F^. The physical meaning of the covariant derivative is 




c -f 



dcosu 



FIG. 2: Schematic representation of the c and c vector defined by the layer normal vector m and 
the director n. It is obvious that c — > — c in the small tilt limit a — > 0. 

as follows (Fig. 2). With the sinusoidal density profile ^(r) = ^Fo exp (i<f)(r)) and the layer 
normal vector m = V0/|V0|, one obtains D^/ = i\V<f>\(m — ncosa)^ = ic\V(f>\^, where a 
is the angle between n and m, and c is the layer normal vector component perpendicular to 
the molecular orientation, which coincides with — c for a — > 0. Because Cy |n- D\l/| 2 = and 
C±\n x D^F| 2 = C±^l\\7(f)\ 2 c 2 , Sm-C phase is equilibrated with a combination of a negative 
Cj_ and stabilizing fourth-order covariant derivatives. 



The scattering function is readily calculated as 



a/ \ Q2) 

% + Diifaii - g ) 4 + (A + 2D\\±qi)(q\\ ~ go) 2 + D±(q 2 ± - g \) 2 ' 

where we set fo = f — Cj_/(4D±) and g 2 j_ = —C±/(2D±). The relations of the present 

generalized model to the other models are shown in Table I. One can obtain the original 

Chen-Lubensky model if D\\ = D ]l± = and D ± ^ 0. With C\\ = J B||/(2^ 2 g 2 ), C ± = 

Bj_/ '(2\E'ogo) an d Al — D\\± = D± = 0, the Chen-Lubensky model is reduced to the Landau- 

de Gennes model (la). 

model type parameter list 

original Chen-Lubensky model [9] Dn = Din = 0, D± / 

model of [42] D\\ = = D ± ^0 

model of [43] (in the vicinity of q\\ = q ) C\\ = ±D\\ql ^ 0, D\\ ± , D± ^ 

Landau-de Gennes model (la) C t = Bj/(2*ggg) (i =||, _L), D\\ = D\\ ± =D ± = 

TABLE I: Relationships between the generalized Chen-Lubensky and other models with the cor- 
responding parameter lists. 

The equilibrium director in the Sm-C phase tilts against the layer normal at the angle, 

(13) 




Next we consider a perturbation of the uniform equilibrium configuration with the gen- 
eralized Chen-Lubensky model. The director perturbation is given by 

5n(r) = n(r) - N, (14) 

where N is the equilibrium director and equals (sin a cos ijj, sin a sin cos a) . The azimuthal 
angle ip is the Goldstone mode in the Sm-C phase. The layer deformation is expressed only 
by u(r) 



\J/(r) = ^ exp 



(15) 



The equilibrium condition (5), equivalent to n || 5F/5n, reads 

5F Ni + 5rii 5F 



5rii N 3 + 5n 3 5n 3 
9 



(i = 1,2). (16a) 



The perturbation expansion requires 



7i 



5F {0) 



5F {1) 



li 



5F 



(i) 



5n,5F {0) 



= (2 = 1,2), 



(17) 



'5n3 5rii 5n 3 N 3 8n 3 

where 7$ = Ni/N${i = 1,2), and 0^ n ^ means the nth order perturbative part of O in terms 

of 5n and u. Equation (17) has two possible solutions. We obtain the Sm-A phase (iV = e 3 ) 

/ ~ \(0) 

if ( 5F/Sn 3 ) 7^ 0, otherwise the Sm-C phase with the director tilt angle a (13). 



B. Effective Sm-A free energy 

The Sm-A case is quite simple: 71 = 72 = and A 3 = 1. Following the same procedure 
as in Sec. II, the effective free energy is calculated as 

2i(C ± + D liq l + D ± q 2 ± ) 



Sni(q) 
F 



2C ± + D'q 2 ± + q 2 



-qMq) 



1,2), 5n 3 (r) = 0, 



(18) 



C\\qf\ + D {l ql + 2D ll± qjq ± + D ± q{ 
C±ql (b'ql + q 2 )-2ql (p^ + D^ 



2C ± + D\\ql + D ± ql 



2C ± + D'q{ + q 



HQ) 



(19) 



where fq[- ■ ■] = / dq[- ■ ■ ]/(2n) 3 and D = 2D± + k±. The twist and bend contributions 
of Frank energy do not appear in the effective free energy because 5rij_ || q ± . The layer 
compression term has no director contribution, because the layer width strain energy should 
be the quadratic form d z u + 5n z , and the higher order term 5n z is neglected. 

Next we compare the previous Landau-de Gennes and the present Chen-Lubensky model 
briefly, to understand the role of the fourth order gradient terms. By setting D\\ = D\\± = 
k\ — 0, Eq.(19) is reduced to the simpler form 

C± + D±q 2 ± 




C\\qj + K(q)ql + 



2C ± + 2D ± q{ + q 

C ± + 2D ± C + D ± qj 
2C ' 



2 



(20) 



2D 1 _q\ + g 2 

This is identical with the effective Landau-de Gennes free energy (7) if D± = 0. The profile 
of k(q) is depicted in Fig. 3. The layer bending modulus only from the density contribution 
can be considered as k(q± — > 00) = 2D_ L /(2D_ L + 1) because the splay vanishes at q± — > 00. 
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1 








K(0) 


=4 







(2d ± +\y 



FIG. 3: Plot of the effective layer undulation modulus h{q^ = 0, q±) in (20). The Lorentzian curve 
with the offset (~ D±) decays with the characteristic wave number ~ \J C±. The fourth-order 
gradient coefficients are typically quite small than the second order coefficients. 

C. Effective free energy for Sm-C phase 

Next we calculate the effective free energy for the Sm-C phase. Using 5F^/8n^ = 
—4 + 2Z)_l<?o ((V-Wf) ~~ 1)} = 0' the perturbed equation of state is given in the 
matrix form 



M5n ± (q) = v(q), 



(21a) 



M = 



v(q) 



2i 



D'q 2 2 + (1 - N 2 ) q 2 + 2N 2 A + k 2 q 2 ±2 -D ' q x q 2 + N t N 2 (q 2 + 2A) - k 2 q ±1 q ±2 
-D'q x q 2 + AqiV 2 {q 2 + 2A) - k 2 q ±1 q ±2 & q 2 + (1 - N 2 ) q 2 + 2N 2 A + k 2 q\ x 

(21b) 

i?7i + dq 2 
E^ 2 - dqi 



where we use the first order terms of the normalization condition Sn^ 
Here we introduced the abbreviations 

2D ll± q 2 N 2 



1 



A = tttt I Cm + 4D ± q 2 + 



E 
Q 

d 




2D\\±q 2 N 2 



Nq, 



Nxq, 



D ll± (N-q) 2 + D ± (Nxqf. 



-7i5ni — 72<5n 2 . 

(22a) 

(22b) 

(22c) 
(22d) 
(22e) 



(21c) 
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The smectic free energy is invariant under a uniform layer displacement, and the solution of 
the equation of state (21a) for q = is N ■ 5nj_ = 5n 3 = 0. This is nothing but a uniform 
rotational Goldstone mode of the c-vector. 

Inverting the matrix M, we obtain the expression for the director deformation 



/ 



5ni(q) 
5n 2 (q) 



2iu(q) 



D 



c 



-N 3 (ED'q 2 + 2dAW 2 ) q 3 + EN 3 N,<; 



2\ 



-En 2 q ±1 q ±3 - N 3 dq ±1 (q 2 + K 2 q 2 ) 
N 3 (ED'qx + 2dAN^j q 3 + EN 3 N 2 q 2 
-Ek 2 q ±2 q ±3 - N 3 dq ±2 (q 2 + k 2 q 2 ) J 



\ 



(23) 



where the determinant of M is given by 



D c = Ni 



D + k 2 ) q 2 q 2 + D'R 2 f + 2D' Aft + 2Ak 2 q z ±3 + 2N'iq z A + q 



,~2 



r2 „2 



(24) 



Now in the Sm-C phase, 8n !ff q L even at the pure undulation (q 3 = 0). So not only the 
splay Frank elastic energy but also the twist and bend terms have contributions to the layer 
bending elastic energy. The splay term favors the state with q _L 5n (we call it the q _L 5n 
state), and both the twist and the bending terms favor the q \\ Sn state. In such a layered 
smectic phase with the spontaneously symmetry breaking layer normal direction and with 
the rotational Goldstone mode of the c-vector, the three Frank elastic terms have anisotropic 
energy contributions, together with the anisotropic coupling constants. The leading order 
terms including D\\, D± and D\\± are (V|w) 2 , (V^u + V -n) 2 and ( V|m) (V5_m + V-n), which 
favor N ± q, N \\ q and the intermediate state respectively. 

Anisotropy due to the c-vector is apparent in the pure undulation case (q 3 = 0). Let i? 
be the angle between c and q ± . The effective free energy is calculated with (23), 



2 1 f = A q{\u(q)\ 2 (K u (x) + K n (q,x)), 

'q± J< i± 



(25a) 



K u {x) 



K, n (q,x) 



D\\Nix 2 + 2D ll± Nlxy + D ± y 2 



N 2 



1 D nl (x)N 2 q 2 ± \ 2 A + D n2 (x)N 2 ± 
N 2 D n3 {x)N 2 ql\ 2 A + D n4 (x)Nl' 



(25b) 

(25c) 
(25d) 
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D nl (x) = -2 (l) ll± x + D ± y^j y (1 + k 2 y) , (26a) 
D n2 (x) = xNl jl + D'y - 4 (b^x + £> ± y) } (1 + k 2 y) , 

-4(b ]l± x + D ±y y (l-x), (26b) 

D n3 (x) = (l + D'y) (l + fay), (26c) 

D n4 (x) = 2 (1 - x) + k 2 N*x + l) , (26d) 

where x = cos 2 -^ and ?/(#) = 1 — ATja;. The effective layer bending elasticity contains the 
two components: k u (x) is the density contribution without the director deformation, and 
k u (q±,x) comes both from the density and director elasticity. The director part has a new 
characteristic length-scale X A = 1/ \J~A, determined by the ratio of the Frank elasticity and 
the combination of the second and fourth order density gradient terms (22a). Taking the 
Landau-de Gennes limit .Dj — > (i =_L, ||_L, ||) and a — > 0, is reduced to the parallel 
penetration length Ay = It is completely different from the Sm-A case where the 

characteristic length is Aj_ = \/(Tl (20). 

We will compare this result with the previous works in Sec. Ill E. 



D. State diagram for Sm-C phase 

In the following, we examine the minimizer d = 8 for F("&) (25a). The angles d and 9 
are spatially uniform by assumption. 

We first consider the one Frank constant case: K 1 — K 2 — K 3 . In this case, the calcula- 
tion and evaluation of the free energy become quite simple with no further approximation. 
One can show that the free energy density f(x) is a convex function throughout < x < 1, 
because a stability of the smectic layer requires the inequality det(-D) = D\\D± — D^ ± > 
(see Appendix B). Thus the behavior of the free energy minimum is determined only by 
/'(0), f (I) and det(-D). There are three possible cases for the first derivative f'(x): 

(i) /'(0) > (so /'(l) > 0) 

The minimum of the free energy is at 9 = 90° 

(ii) /» < and /'(l) > 

The favored angle 9 exists between 0° and 90°. 

(iii) /» < and /'(l) < 0. 
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The free energy minimum is at 6 = 0°, and the director tends to be aligned to the layer 
bending direction. 



(a) a = 10°, q ± =2x lO^m" 1 (b) a = 10°, q± = 2 x lOW -1 




FIG. 4: State diagrams with D±, D\\± and D\\ in the one Frank constant approximation. We use 
q = 2 x 10 7 cm -2 and Cy = 4 x 10 12 cm _1 Ajj~ 2 ^ as a set of typical values [18, 44]. The planes 
/'(0) = 0, /'(l) = and det(£») = - D^ ± = are denoted by the plane (1), (2) and (3) 

respectively. 

The fixed parameters are set to be typical values and used in followings if not specified: 
q_i_ = 2 x 10 7 cm _1 and C\\ = X^ 2 = 4 x 10 12 cm~ 2 [18, 44]. We limit the possible parameter 
range of |.Dj | i =\\, ||_L, J_| to < Di < 1, because of the stability of the free energy, 
the relation D± ~ d 2 /(\±N±) 2 holds, and the feasibility of the Landau-de Gennes model 
(Di = 0) in the Sm-A phase. In fact, one of the scattering experiments shows that D\\±, 
D\\ ~ and D±(~ 10 -4 ) 1 [15], where only the trivial 9 = 90° state is expected to 
be observed as the previous work implies [30, 35]. However, depending on the material, 

14 



one might be able to obtain a wide range of the parameter sets within < £>i < 1. The 
three-dimensional state diagrams for D\\, D\\±_ and D± are given in Fig.4. State boundaries 
are determined by the conditions (1) /'(0) = 0, (2) /'(l) = 0, and (3) det(L>) = 0. 

We again note that the smectic layer elasticity consists only of the density and the Frank 
elastic energies and the coupling energy between the layer normal and the director. All the 
following results will be explained in terms of combinations of the three elastic contributions. 

The 9 = 90° state dominates for small a and q± < A^(~ 10 7 cm~ 1 ) (Fig.4(a)). When 
Q± ^ A^ 1 , on the other hand, the Frank elasticity less contributes to k(q), and anisotropic 
density terms dominate (Fig.4(b)), being enhanced for larger a (Fig. 4(c)). Thus various 
angle 9 can appear. In experiment, the largest tilt angle a that has been reported to our 
knowledge is 37° [45]. So the possible stable 9 state might be limited to 0° < 9 < 90°, as we 
see from the diagrams for a < 40°. Nevertheless, if the a = 50° molecular tilt is achieved, 
there exists the stable 9 = 0° state (Fig.4 (d)). 

In each state diagram, the behavior of 9 is quite sensitive to D±_ compared with D\\± 
and D\\. This is because the elasticity of the vector c (~ — c) is dominated by the elastic 
terms including D±, as we discussed in Sec. II A. The ratio between C± and D± is uniquely 
determined by the tilt angle a (13), so the value of D± controls the preference of q ± \\ N. 
Thus the stable state changes from 9 = 90° to 0° < 9 < 90° and 9 = 0° as D ± grows. On 
the other hand in Fig.4(d), the 9 = 0° state transits to the 0° < 9 < 90° state with the 
increment of D\\. 

The state diagram for the reduced Frank elastic moduli k\ and k 2 is also calculated based 
on (25a). We use D± = 0.5, D\\± = 0.01 and Dn = 0.1. An experimental result tells that 
is smaller than Dj_, D\\ [15]. Dj_ = 0.5 is about on the state boundary between 9 = 90° 
or less, so that we would obtain the sensitive state behavior. This parameter set will be 
used in the following if not specified. 

Effect of the Frank elastic constants is quite sensitive to a. For a < 40° (Fig.5(a)), the 
9 = 90° state dominates except for low k 2 , while 9 < 90° domain shrinks for smaller and 
larger k\. The 9 = 90° state is achieved by a combination of the locking term (which is 
the major contribution for small a), and the twist Frank term \N ■ (q ± x 5ri) | 2 favoring 
Q± II N or q ± \\ 5n. The 9 < 90° domain is suppressed by the two mechanisms. One for 
larger k\ is a combination of the locking term and the splay Frank term (see Fig. 7 and the 
detailed calculation is found in Appendix C). The other for smaller k\ is a combination of 
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(a) q = 40° 



(b) a = 50° 



(c) a = 50°, D± = 0.2 




-L 
j 


r 




9=90° 


1 <0 <90 - 

t 



10 



2 4^6 8 10 




FIG. 5: State diagrams in terms of hi and k 2 with different tilt angles, (a) The 0° < 6 < 90° 
domain gradually disappears as a — > 0. For a = 40°, the 0° < 9 < 90° domain grows around 
h\ ~ 0. (b) For a = 50°, a stable = 0° domain emerges at k\ ~ and higher K2- The preferred 
angle changes discontinuously between the 9 = 0° and = 90° states. The lower D± diagram is 
plotted in (c). 



(a) a = 20° (solid line) (b) a = 50°, ki = 0, R 2 = 6 

a = 40° (dashed line) (K1/K3 = 1, ^2/^3 = 7) 

1.2 




(deg.) 



(deg.) 



FIG. 6: Plot of the free energy /(??) as a function of the angle between the c- vector and q^. (a) 
The tilt angle dependence of the free energy /($). (b) Discontinuous transition is explained with 
the two minima at 9 = 0° and 90° as it should be. 



the twist and bend Frank elastic energies, which favor q ± \\ SN. For a > 50° (Fig.5(b)), 
the locking term less contributes and the anisotropic elasticity with D± is enhanced more. 
Thus both of the two minima q ± || N and q ± || 5n are stable for higher twist k 2 . They are 
incompatible because of the normalization condition N ■ Sn = 0, so the transition between 
9 = 0° and 9 = 90° is discontinuous. This anomalous state behavior disappears for weak 
anisotropy (Fig. 5(c)). 

We also plot the free energy /($) in Fig. 6. The normalization is taken with /(■$ = 0°). 
The minimizer 9 decreases as a grows, but the free energy barrier is still small for a = 40°, 
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(a) = 0' 



(b) = 90' 




FIG. 7: Layer undulation and the director configuration assuming the nonzero tilt angle a. The 
splay Frank elastic energy costs at d = 0° more than at $ = 90°, because of the director component 
projected to the zx-plane. 

being a few percent of the absolute value Fig.6(a). So the equilibration to the free energy 
minimum might be disturbed by the internal thermal noise and an external force. However 
the 0° < 9 < 90° state is more stable at higher a. The discontinuous 9 transition in Fig. 5(c) 
is certainly due to the coexistence of the two local minima (Fig. 6(b)). 

E. Discussion 

We next compare our calculation with the previous works. Hatwalne and Lubensky [30] 
derived the elastic free energy for the Sm-C phase in a more elegant way making use of the 
symmetry and covariance of the model free energy. The resulting representation in the case 



We can obtain the layer compression coefficient from (11a), (23), and the normalization 
N • 5n = 0, by setting q 1 = q 2 = and q 3 — > 0, 



We note that the layer compression elasticity stems from the director contribution R\ as 
well as the density contribution in the Sm-C phase, because a change of the layer width is 
associated with the change of the tilt angle. 

As for the layer bending elasticity, the effective Chen-Lubensky model (25a) cannot be 
written in the Hatwalne-Lubensky g-dyadic form (27), because in general the anisotropic 



N = (N u 0,N 3 ) is 




(27) 




(28) 
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layer bending elasticity should be expressed as 

F = lj q ( B <ll + <limQ±))Wq)\ 2 . (29) 
Still we can give some representations for Ku and K 22 in the simple form, 

K n = K{d = 0°) = || [2D\\Ni + iV 3 2 (l + R.Nl) } , (30a) 

K 22 = K(0 = 90°) = 2 f 3 ^ {kl + (30b) 

NI (2D ± + «x + lj 

i^i 2 is calculated with the parabolic approximation, 



#(0 = 0°) + = 90°) 1 


^<9K 


dK 


2 4 


^ dx 


i?=0° ^ 



(31) 



The three elastic constants K n , Ku, and i^ 2 2 contain both the density and the director 
contributions as pointed out in the previous research [30]. 

In this way, we bridge the models of smectic liquid crystal at different coarse-graining 
levels. It leads one to a quantitative discussion with the phenomenological macroscopic 
model [30], based on the experimental determination of the Chen-Lubensky model [15]. 

We conclude this section comparing our results for Sm-C liquid crystals with other layer 
forming materials. Diblock copolymers, chemically connected two different homopolymers, 
do not possess the definite molecular orientational degree of freedom especially in the weak 
segregation regime [29]. Thus anisotropic elasticity cannot occur. On the other hand, in 
surfactant systems, the linear amphiphilic molecule is basically normal to the surfactant 
layer [28, 46]. This is similar to the Sm-A phase of liquid crystals rather than the Sm-C 
phase. Thus amphiphilic system again would not have an in-plane anisotropic coupling. 
However in some cases the surfactant molecules may be aligned at a nonzero tilt angle 
against the layer normal [47-49]. Such systems might have an anisotropic layer elasticity as 
in the present Sm-C liquid crystal case. 



IV. HYDRODYNAMICS OF SMECTIC-C LAYERS 

In this section, we consider the transverse wave in smectic liquid crystals which oscillates 
in the layer normal direction, and travels perpendicularly to the layer normal. The wave 
frequency is much lower than the inverse of the molecular time scale (~ 10 8 Hz), so that the 
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permeation effect which changes the layer width is negligible [1]. Also we eliminate some 
rapidly relaxed degrees of freedom, namely the density p, the local energy e, the velocity 
along the layer in-plane directions, the pressure P, and the c-vector in the time scale ~ 10~ 4 s. 

We derive the equation of motion of u(r) from the time evolution of the velocity field 
v(r). According to the unified hydrodynamic description [1, 31], a set of the hydrodynamic 
equations at a constant temperature with no external field is 

V-v(r) = 0, (32a) 

^-..<r) = 0, (32b) 

-B z fadiVp(r) = -vh z {r), (32c) 

Po ~W^ = ~Ju ~ Bzz P d P hz ( r "> +Vz^sd y d s v l3 {r), (32d) 

where h z is the torque acting on the c-vector and B a ^ is the coupling constant between the 
torque and the velocity. The first equation is the incompressibility condition originated from 
the equation of continuity. The second and third equation describes the relaxation of the 
layer permeation and the rotational c-vector mode respectively. The last equation is derived 
from the momentum conservation essentially equivalent to the Navier-Stokes equation. From 
the above equations, we derive the equation of motion as 

d^uiv) du(v') 5F 

po^t- = md^— - -. (33) 

In the smectic-C phase, the model free energy has the static elastic coupling of the layer 
displacement field u(r) and c-vector. Here exists the dynamical viscous coupling of the 
velocity and the molecular orientation in the first term of (33). Anisotropic viscosity tensor 
[/,py (f3, 7=1, 2) depends on the c-vector, as 

= /i 5/3 7 + /iin /3 (r)n 7 (r), (34) 

where liq and ji\ are constants. The isotropic and rotational viscosities are reported to 
be of the same order (~ 1 Poise) in the experimental studies [45, 50]. We then add an 
oscillatory source term s{r,t). We consider the oscillating wall and cylindrical sources 
s(r,t) = Socos(ut)5(x) and sq cos(ut)5(x)5(y) respectively, to examine the anisotropy in 
the Sm-C phase (Fig. 8). This vibrating force can be provided by an oscillatory object, for 
example, an electric field [51]. The linearized dynamic equation in Fourier space is 
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(a) (b) 




FIG. 8: Illustration of the oscillatory sources. We consider the two geometries: (a) the one- 
dimensional wall and (b) the two-dimensional cylinder at the center point of the layer plane. 



Po- 



dt 2 



d 

Oo 5 Pl + ^NpN^q^— + (Bqj + K{$)q\) 



u(q,t) + s(q,t), (35) 



dt 

where s(q,t) = 4ir 2 s cos(ujt)5^(q ± ) and 27rs cos(ujt)5(q z ) is the one-dimensional wall and 
two-dimensional cylindrical source respectively. The thermodynamic force SF/ 5u* is written 
in the linear form L(q)u(q, t) = (^Bq 2 + K("&, q±)q*±j u(q, t). 

We next derive the dimensionless dynamic equation. The characteristic parameter set is 
listed in Table II. The reduced dynamic equation is 

TABLE II: Parameters and their typical scales in the dynamic equation (35) [1, 45, 50, 52]. 

parameter description typical scale 

po averaged molecular density ~ lg/cm 3 

u(r) layer displacement field ~ 10~ 7 — 10~ 8 cm 

Ho isotropic viscosity ~ 1 Poise 

fii rotational viscosity ~ 1 Poise 

qi undulation wave number ~ 10 _5 cm 

u) undulation frequency ~ 10 4 Hz 

B = Cj^Q^Q (i =_L, ||) layer compression elastic constant ~ 10 6 dyn/cm 2 

(second order density elastic constant) 

A^o^o (* = -L) -L|| 5 II) layer bending elastic constant ~ 10 -6 dyn 

(fourth order density elastic constant) 

Ki (i = 1, 2, 3) Frank elastic modulus ~ 10 -6 dyn 
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M 4 d 2 u(q,t) i 2 r] c , „ 2 2 x du(q,t) 



Bqi + K(#)q* u(q,t) + ^s(q,t), (36) 



>-3 



where we introduced the dimensionless valuables t = rt, lot = Co, uq = u, q£ = q iy 
Vi — VcVi (i — 0? 1)) B = B B and K(-&) = K 3 K($). All the dimensionless constants are of 
order ~ I. We set the characteristic scales according to Table II: r = lCT 4 s, £ = 10~ 5 cm, 
7] c = 1 Poise and B = 10 6 dyn/cm 2 . In this parameter set, the scale of each term in (36) is 
p ^ 4 /K 3 r 2 ~ 1CT 6 and £, 2 r} c /K 3 T ~ B /K 3 q 2 ~ 1. Thus the inertia term can be neglected 
and the final form of the dimensionless dynamic equation is 

-du^t) + ~ L (q)u(q,t) = ~s(q,t), (37) 
ot 

where the effective viscosity ( = q 2 (fjo + fji cos 2 $) , and the normalized external force s 
is given by s = s (q) cos(Cot) . The reduced amplitude s is ^ 7l2s o ( lo^^( ( ly)^( ( lz)/K 3 and 
2nsoqo£ 4 '5( y q z ) / K 3 in the oscillatory wall and the vibrating cylinder case respectively. The 
time evolution is 

u(q, t) = —fz — r- ( — ^ cos Cot + Co sin Cot + A exp At J , (38) 

C(A 2 + d> 2 ) V ) 

where A = —L/C, and the initial condition is set to be u(q, 0) = 0. 

Visualization of the wave propagation is depicted in Fig.9 and Fig. 10. Unphysical infinite 
wavelength mode is excluded. We use the parameter set: C\\ — 4 x 10 12 cm -2 , D± = D\\± = 
D\\ = 0.1, Hi = k 2 = 0, uo = 5 x 10 4 Hz, and a = 30°. The reference c-vector is tilted against 
the x-axis at 40°. The same parameter set is used in the following if not specified. The layer 
displacement is shown with the brightness, and the director change is indicated with the 
unit arrow. The dumped director degree of freedom is governed by the layer displacement 
field through (23). While the layer displacement scale is ~ 10~ 7 cm _1 , the director rotation 
angle is about 10~ 3 — 10~ 4 rad, which can be observed in experiment [53, 54]. 

In Fig.9, the planar oscillation orientates the director with the layer displacement. 
Anisotropic wave propagation is obvious under a cylindrical source (Fig. 10). The favored 
angle is 9 — 90°, as in the experimental study [35] and the free energy analysis (Fig. 4). The 
wave velocity is faster in the rigid direction c || q ± , than in the soft orientation c _L q ± . A 

simple dimensional analysis of (37) with 5 = gives the ratio 

i 

v\\/v ± = (X(0°)/X(90°)) 5 {f!o/ (fj + 771)}* • (39) 
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We conclude this section with the remark that both the director alignment and the 
anisotropic wave propagation are entirely controlled by the static and dynamic coupling 
between the layer displacement and the director (39), and could be observed within the 
experimental resolution. 
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FIG. 9: Snapshots of the layer oscillation under a planar source (x = 0). The brightness indicates 
the local layer displacement and the unit arrow represents the deviation of the director from the 
ground state. We set C\\ = 4 x 10 12 cm~ 2 , D ± = D\\ ± = D\\ = 0.1, k ± = k 2 = 0, u = 5 x 10 4 Hz, 
and a = 30°. 



V. CONCLUSION 

In this paper, we analytically derived the effective generalized Chen-Lubensky model by 
adiabatic elimination of the director relaxation. For Sm-A phase, the director unlocking 
from the layer normal is confirmed at an undulation wavelength shorter than the director 
coherent length. The effective layer bending elastic modulus is written as a function of the 
unlocking angle a, and decays with the increase of a. This agrees with the argument in the 
previous work [25]. In Sm-C phase, on the other hand, not only the director unlocking but 
anisotropic elasticity arises from the c-vector. After the detailed study, it turned out that 
the preferred angle 6 between the layer bending orientation and the c-vector varies from 0° 
to 90° depending on the layer elastic constants, the Frank constants, the tilt angle, and the 
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FIG. 10: Anisotropic wave propagation is shown in the presence of an isotropic oscillatory source 
(x = y = 0). The wave front is elliptic with the long axis parallel to the reference director c. 



undulation wavelength. Transitions between the different 6 states are not only continuous 
but also discontinuous. Then the new characteristic length is found, which plays an 
important role on the elasticity of Sm-C phase. The model parameters of the macroscopic 
free energy [30] are determined from the more microscopic Chen-Lubensky model. It allows 
a quantitative argument based on the model [30] with the experimental determination of the 
model parameters [15]. We also compared smectic liquid crystals with other layer forming 
materials. 

Next we discussed the hydrodynamics of Sm-C layers using the effective elastic energy 
derived above. We derive the Sm-C hydrodynamics with the director deformation elimi- 
nated. Anisotropic wave propagation from the isotropic oscillatory source and the director 
orientation in the simple undulation along one direction are confirmed. These effect could be 
observed in experiment [51, 53, 54]. By converting the mechanical undulation to the optical 
information through the director orientation, here also arises a new possible application of 
liquid crystals for a sensitive mechanical sensor. 
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APPENDIX A: EFFECTIVE LAYER BENDING ELASTIC MODULUS AS A 
FUNCTION OF m-n UNLOCKING ANGLE 



FIG. 11: Schematic representation of the smectic layer with a sinusoidal undulation. Here exists 
the small unlocking angle a between the layer normal m and the director n, which accounts for 
the weak bending elastic modulus. 

Here we give a better representation of the effective layer bending elastic modulus (9), 
as a function of the root mean square (RMS) of the angle a between the layer normal and 
the director. This leads us to a more quantitative understanding of the previous numerical 
result on the mean curvature of smectic layers in the TGB phase [25] . 

We suppose a single Fourier mode of the layer undulation u(r) = u cos(qx) (Fig. 11). 
The wave-number vector is along the x-axis. With the help of (6), m and n are 




1 



qu sm(qx), 



(Ala) 
(Alb) 



n. 



'X 



l + (Ag) 2 " 
quo sm(qx), 



m 



'X 



m y = 0, 



(Ale) 



m z = 1. 



(Aid) 




qu sm(qx). 



(A2) 
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The RMS a is readily calculated 



1 u X 2 q 3 1 2 3 
ot = — ^- — ~ —=u X q , 



V2l + (\q) 2 ~ ' (A3) 

where the smooth undulation condition q <C A -1 is assumed in the last line. We note that 

the tilt angle is roughly proportional to the cube of the wave-number, and the characteristic 
length is a combination of the penetration length and the undulation amplitude (A 2 m ) 1 ^ 3 . 
Thus the effective layer bending rigidity is now written as a function of a 

1 



k, (a) 



2/3 ' 



(A4) 



1 + (V2a\/u ) 

The same method could be applied for the Sm-A phase also in the Chen-Lubensky model. 



APPENDIX B: DERIVATION OF L> ± L>| - D 2 ± > AND ONE CONSTANT AP- 
PROXIMATION 

Here we prove the relation D±D\\ — D 2 ± > to guarantee the convexity of the effective 
free energy as a function of x — cos 2 i9 in the one Frank constant case (B4). This relation 
actually holds for an arbitrary set of Frank elastic moduli, because the inequality is based on 
the stability of the smectic layer. With the density profile ^ = exp(qr-r), the g-dependent 
part of the Chen-Lubensky Hamiltonian (11a) is proportional to 

DuQj + 2D u QlQ 2 ± + D L Qi 



where the effective momenta Q 2 
characteristic wave number q 2 „ = 



(Bl) 



(Q\\ ~ <?o) 2 + <?o|| and Q\ = q\ - q 2 ± , the 



C\\D ± -C ± D» ± )/ 



2 DiDii - 2Du 



C\\D\u-c ± dA/ 



2[D\D\\ -2Dl 



and q 2 ± = 

and the projected wave-number vector = 
(n ■ q)n, q L = —n x (n x q) are utilized. Positive definiteness of (Bl) requires 

D ± D\\ - Dj ± > 0. (B2) 

This inequality helps one to prove a convexity of the free energy at the one constant ap- 
proximation (Ki — K 2 — K 3 ). In this case, the free energy is simplified as 

_) x 2 + D ± (1- Nix) + 2D\\±N]_x 

(2N 2 ± + k\) (-2NiD 2 ± x 2 + D ± (1- Nix) + 2D\\ ± Nlx} - 2N 2 NlD ± x 



q\\u{q) 



2iVf [ D\\D± - Dl ) x 2 



2k 2 D\ (1 — N 2 x) + 4D ± N 2 (l-x) + 2N 2 



(B3) 
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and the second derivative is 

-J\x) = D {l D ± -Df u 



N 2 D A 



2N 2 ± + k\ 



ANiq\_\u{q)r 



{ (l + 2D^j (2N\ + k\ ) - 2N 2 ± D ± (k 2 ± + 2) 



x { (l + 2D^j (2Nl + k\) D l]± - (k\ + 2) D ± Y , 
where the dimensionless wave number is given by k±(D\\±, D±) = \aQ±- 



(B4) 



APPENDIX C: SPLAY FRANK ENERGY AS A FUNCTION OF 

We calculate the splay energy for a small tilt angle a and the angle between c and 
q ± . Let the rr-axis be parallel to q ± . The director is n = (sin a cos sin a sin cos a). 
Assuming the smooth layer undulation q± A -1 , the director deformation perfectly follows 
the layer normal vector. Undulated director configuration is obtained by operating the 
rotation matrix about the y-axis R y 



^ cos qx — sin qx ^ 



Ry 



1 

^ sin qx cos qx j 

to the reference director n. Thus the spatially averaged splay energy is expressed as 

a 2 

{V • (Ryn)} 2 ) = — (sin 2 a cos 2 d + cos 2 a) . 
2 



(Cl) 



(C2) 



Stability of the 9 = 90° state grows as the splay Frank elasticity dominates. However this 
effect of the splay term is not very strong, due to the factor sin 2 a. 



[1] P. G. de Gennes and J. Prost, The Physics of Liquid Crystals (Clarendon, Oxford, 1994). 

[2] P. M. Chainkin and T. C. Lubensky, Principles of Condensed Matter Physics (Cambridge 

University Press, Cambridge, England, 1995). 

[3] P. H. Keyes, H. T. Weston, and W. B. Daniels, Phys. Rev. Lett. 31, 628 (1973). 

[4] A. Aharony, R. J. Birgeneau, J. D. Brock, and J. D. Litster, Phys. Rev. Lett. 57, 1012 (1986). 

[5] R. Shashidhar, B. R. Ratna, and S. K. Prasad, Phys. Rev. Lett. 53, 2141 (1984). 

[6] Y. Drossinos and D. Ronis, Phys. Rev. A 33, 589 (1986). 

26 



[7] P. G. de Gennes, Mol. Cryst. Liq. Cryst. 21, 49 (1973). 

[8] K. C. Chu and W. L. McMillan, Phys. Rev. A 15, 1181 (1977). 

[9] J. Chen and T. C. Lubensky, Phys. Rev. A 14, 1202 (1976). 
[10] L. Benguigui, J. Phys. (Paris). Colloq. 40, C3 (1979). 
[11] C. C. Huang and S. C. Lien, Phys. Rev. Lett. 47, 1917 (1981). 
[12] G. Grinstein and J. Toner, Phys. Rev. Lett. 26, 2386 (1983). 
[13] S. Witanachchi, J. Huang, and J. T. Ho, Phys. Rev. Lett. 50, 594 (1983). 
[14] C. R. Safinya, L. J. Martinez-Miranda, M. Kaplan, J. D. Litster, and R. J. Birgeneau, Phys. 

Rev. Lett. 50, 56 (1983). 
[15] L. J. Martinez-Miranda, A. R. Kortan, and R. J. Birgeneau, Phys. Rev. Lett. 56, 2264 (1986). 
[16] H. S. Kitzerow and C. Bahr, eds., Chirality in Liquid Crystals (Springer- Verlag, New York, 
2002). 

[17] P. G. de Gennes, Solid State Commun. 10, 753 (1972). 

[18] S. R. Renn and T. C. Lubensky, Phys. Rev. A 38, 2132 (1988). 

[19] J. W. Goodby, M. A. Waugh, S. M. Stein, E. Chin, R. Pindak, and J. S. Patel, Nature 337, 
449 (1989). 

[20] T. C. Lubensky and S. R. Renn, Phys. Rev. A 41, 4392 (1990). 

[21] M. Ismaili, A. Anakkar, G. Joly, and N. Isaert, Phys. Rev. E 61, 519 (2000). 

[22] E. Grelet, B. Pansu, M. Li, and H. T. Nguyen, Phys. Rev. Lett. 86, 3791 (2001). 

[23] B. A. DiDonna and R. D. Kamien, Phys. Rev. Lett. 89, 215504 (2002). 

[24] J. Yamamoto, I. Nishiyama, M. Inoue, and H. Yokoyama, Nature 437, 525 (2005). 

[25] H. Ogawa and N. Uchida, Phys. Rev. E 73, 060701 (R) (2006). 

[26] E. Grelet, B. Pansu, and H. T. Nguyen, Phys. Rev. E 64, 010703(R) (2001). 

[27] T. Ohta and K. Kawasaki, Macromol. 19, 2621 (1986). 

[28] G. Gompper and S. Klein, J. Phys. II (France) 2, 1725 (1992). 

[29] L. Leibler, Macromol. 13, 1602 (1980). 

[30] Y. Hatwalne and T. C. Lybensky, Phys. Rev. E 52, 6240 (1995). 
[31] P. C. Martin, O. Parodi, and P. S. Pershan, Phys. Rev. A 6, 2401 (1972). 
[32] A. Buka and L. Kramer, eds., Pattern Formation in Liquid Crystals (Springer New York, 
1996). 

[33] A. N. Pargellis, P. Finn, J. W. Goodby, P. Panizza, B. Yurke, and P. E. Cladis, Phys. Rev. A 

27 



46, 7765 (1992). 

[34] T. Carlsson, F. M. Leslie, and N. A. Clark, Phys. Rev. E 51, 4509 (1995). 

[35] D. Johnson and A. Saupe, Phys. Rev. A 15, 2079 (1977). 

[36] P. E. Cladis, Y. Couder, and H. R. Brand, Phys. Rev. Lett. 55, 2945 (1985). 

[37] C. Chevallard, T. Fisch, and J. M. Gilli, J. Phys. II (France) 7, 1261 (1997). 

[38] N. Clark, Appl. Phys. Lett. 35, 688 (1979). 

[39] N. Clark and R. B. Meyer, Appl. Phys. Lett. 22, 493 (1973). 

[40] S. V. Yablonskii, K. Nakano, M. Ozaki, and K. Yoshino, JETP Lett. 77, 140 (2003). 

[41] S. Uto, E. Tazoh, M. Ozaki, and K. Yoshino, J. Appl. Phys. 82, 2791 (1997). 

[42] I. Luk'yanchuk, Phys. Rev. E 57, 574 (1998). 

[43] A. Kundagrami and T. C. Lubensky, Phys. Rev. E 68, 060703 (2003). 

[44] L. Navailles, R. Pindak, P. Barois, and H. T. Nguyen, Phys. Rev. Lett. 74, 5224 (1995). 

[45] A. Meiboom and R. C. Hewitt, Phys. Rev. Lett. 34, 1146 (1975). 

[46] K. Chen, C. Jayaprakash, R. Pandit, and W. Wenzel, Phys. Rev. Lett. 65, 2736 (1990). 

[47] S. A. Safran, M. O. Robbins, and S. Garoff, Phys. Rev. A 33, 2186 (1986). 

[48] Z. Wang and C. Gong, Phys. Rev. B 54, 17067 (1996). 

[49] H. Kuhn and H. Rehage, Phys. Chem. Chem. Phys. 2, 1023 (2000). 

[50] B. Tamamushi, Rheol. Acta 13, 247 (1974). 

[51] P. Pieranski, L. Beliard, J. -P. Tournellec, X. Leoncini, C. Furtlehner, H. Dumoulin, E. Riou, 

B. Jouvin, J.-P. Fenerol, P. Palaric, et al., Physica A 194, 364 (1993). 

[52] H.-Y. Chen and D. Jasnow, Phys. Rev. E 61, 493 (2000). 

[53] Y. Galerne, Phys. Rev. A 24, 2284 (1981). 

[54] L.-Y. Kuo, K.-L. Lee, and J.- J. Wu, Jpn. J. Appl. Phys. 45, 8775 (2006). 



28 



